Non-exponential relaxation in disordered complex systems 
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We have analytically obtained the non-exponential relaxation function for disordered complex 
systems applying the multi-level jumping formalism to the fluctuation quantity which makes diffusive 
motion stochastically in the disordered complex space. It is shown that the relaxation function of 
disordered complex systems decays obey to stretched exponential law. 
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Non-exponential relaxation in disordered complex sys- 
tems is the object of active research due to its impli- 
cations in the technology and in several fields of scien- 
tific knowledge. The theoretical and experimental stud- 
ies (for a recent review, see Ref. ^ij) show to us that the 
relaxation behavior in the disordered complex systems 
such in glasses an supercooled liquids, liquid crystal poly- 
mers, dielectrics, magnetic systems, amorphous semicon- 
ductors, pinned density wave, protein dynamics, protein 
folding, and population dynamics among others deviate 
considerably from the exponential Debye pattern Q 



(1) 



and it obey to the stretched exponential relaxation func- 
tion experimentally 



$(i) = *oexp[-t/T]' 



< a < 1 



(2) 



often referred to as the Kohlrausch- William- Watts 
(KWW) function 

The important task in the relaxation in disordered 
complex systems is to extend the Debye theory of relax- 
ation of polar molecules to fractional dynamics, so that 
empirical decay functions, e.g., the stretched exponential 
of KWW, may be justified. In order to explain the origin 
of the non-exponential relaxation in the disordered com- 
plex systems several models, for example, continuous- 
time-random- walk models (CTRW) M S i, "o", To", 111 , 
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models [11 El E3l> defect models 
2^. j um p relaxation models [U , ran- 
dom barrier models |22 . l2a . uM . l25l | , hierarchical models 
^6, 27, 28, 29] have been used and many generalizations 
El El El SEE El Hi E3 of the Debye theory 
have been suggested up to now. 

We remark that one of the most plausible model in 
among others is definitely phase-space model in which 
the possible origin of non-exponential relaxation is ex- 
plained based on the idea of the energy landscape and 
nontrivial energy barrier so that these models provide a 
direct link between the phase-space dynamics and slow 
relaxation. Indeed, several experiments and computer 



simulations have been done which support the explana- 
tion of these relaxation phenomena in the framework of 
energy landscape paradigm as the result of activated dif- 
fusion through a rough energy landscape of valleys and 
peaks 41] . Therefore, it is suggested that such en- 
ergetic disorder which produces obstacles or traps which 
delay the motion of the particle and introduce memory 
effects into the motion 18]. Hence, in such scenario a 
process can be characterized by the temporally nonlocal 
behavior H]. 

Considering stochastic framework, in a relaxation pro- 
cess, the fluctuation variable x which may represent the 
physical quantity such as dipole, jumps at random from 
one value to another with equal probability and it takes 
stochastic values as xi, X2, X3,...,xn with time. Such a 
stochastic process can be regarded as multi-level jump- 
ing process. In this study, using multi-level jumping for- 
malism we shall obtain relaxation function of disordered 
complex system in which Brownian quantity diffusive as 
stochastically. The relaxation function for a stochastic 
process is given simply 



$ (i) = ^' (i = cx)) - ^' (t) 



(3) 



where ^ {t) is the response function. The response of 
the any system can be written in terms of the correlation 
function as 



{x{0)x{t))) 



(4) 



'Electronic address: "ekre m. aydinerQdeu.edu.tr| 



where k is the Boltzmann constant, T is the temperature, 
and (x^) is the mean-square average of the x. On the 
other hand, {x (0) x (t)) is the correlation function of the 
fluctuation variable x. The correlation function can be 
given in discrete form follows 

{x {0)xit)) =J2J2p (^o) ixo\x,t)x (5) 



where the quantities xq and x are the values of stochastic 
variable x at times and t, respectively. Having started 
at the initial state xg with the statistical weight p(xo), 
P {xq \ X, t) measures the probability of propagation from 
Xq to X in time t, which is called as conditional probabil- 
ity. The correlation function measures the decay of the 
fluctuation variable of the physical quantity in the sys- 
tem, and it is determined due to conditional probability 
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which depends on the nature of the stochastic system i.e, 
the system is whether ordered or disordered complex. 

Eqs. (3-5) clearly indicate that relaxation function can 
be obtained from fluctuation-dissipation theory in the 
framework linear-response theory if correlation function 
belong to the relaxation process is known. 

The conditional probability in Eq. (5) for disordered 
complex space may be written formally |43l | , but sugges- 
tive notation, as 



cess) with N values [46 



P{xo \x,t) = d^Q 



1 



(t-t'T 



r(a) JO 

Y,W{x' \x)P{xo\x',t')dt' (6) 



where a is the fractional order which plays the role of a 
dynamical exponent, dxo is the initial condition at t ^ 0, 
and W (x \ x') is the jump rate for a particle from x to 
x' . The Eq. (6) is known the fractal time master equation 
which defines non-Markovian stochastic process with a 
memory, which is a special case of the CTRW equation 

lillii- 

In the operator formalism [il . lisl l46l | the conditional 
probability P (xq \ x,t) is simply given with the matrix 
element of operator P{t) as 



P{xo\x,t) = {x\P{t)\xo). (7) 

where the conditional probability satisfy 

P {xo I X, 0) EE (a;| P (0) |a;o) = 5 (x - xq) (8) 

for t = Q. Similarly, W is also the jump operator which 
is given by 



W{x' I x) = {x\W\x') . 



(9) 



Matrix representation allows to us to write the physical 
quantities as an operator. Hence, the correlation function 
in Eq. (5) can be transformed to the matrix representa- 
tion 



{x (0) X (t)) = (^^oI ^ \xo) {x\ P (t) \xo) {x\ X \x) 

X 

(10) 

where X is the matrix representation of the fluctuation 
quantity x, P{t) is the matrix form of the conditional 
probability, and |...) indicates vector space of stochastic 
states. 

The fact that x takes N valued stochastic variable in 
which the stochastic states have been represented in N 
dimensional vector space. Accordingly, the fluctuation 
variable x and conditional probability P (xq \ x,t) may be 
given by the matrix representation associate a stochastic 
state |x) (x = 1, 2, 3, A'^ for multi-level jumping pro- 
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(11) 



The stochastic states |x) (and similarly |x')) is taken to 
form an orthonormal set which provide closure property 



as 



E 



X) (X 



(12) 



for discrete variables, and assign to it an apriori occupa- 
tion probability p{x) — jj . 

Hence, the fluctuation variable x can be also given 
by matrix form, the eigenvalues of which correspond to 
stochastic variables xi, X2, X3,...,XAr, respectively. The 
matrix form of the operator X of the fluctuation variable 
X is represented as 



X = 



/ xi 

X2 
X3 



V 



\ 







XN J 



(13) 



Now, after these deflnitions, the correlation function 
Eq. (6) can be carried out using this formalism, however, 
we need the operator form of the conditional probability 
in order that calculation of it. 

The time derivative of Eq. (5) is written as 



x,t) 

E 



1 d 
r (a) dt 



/\ a — 1 



W {x' \x)P{xQ\x',t')dt' (14) 



It is seen that Eq. (14) contains a convolution inte- 
gral with a slowly decaying power-law Kernel AI (t) = 
j.a-1 (^a), ensures the non-Markovian nature of the sub- 
diffusion process defined by the fractional diffusion pro- 
cess. This convolution integral is defined as an operator 
which is known as fractional Riemann-Liouville integro- 
differential operator |4Z|; 



oDl-"P{xo \x,t)^ 



1 d 
r (a) dt 



dt 



, P{xo\x',t') 
(t-t')''" 



(15) 



with < a < 1. 

Eq. (14) can be reduced to simply form as down 



,<) =0 A'""E^(^'l^)^(^ol^''*') 



(16) 
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using Eq. (15). On the other hand, using operator formal- 
ism in Eqs. (7) and (9), Eq. (16) can be rewritten com- 
pletely matrix form 

I {x\ P (t) \xo) =0 Dl-- J2 ^ 1^') P W l^o) • 

(17) 

Using closure property in Eq. (12), Eq. (17) is reduced 



to 



{x\ \xo) =0 Dl-^ {x\ WP (t) \xo) . (18) 

This equation can be simplify as 



|p(t) =0 Dl-"WP(t). 



(19) 



Eq. (19) is known fractional relaxation equation li^. 
The solution of it may be presented in terms of Mittag- 
Leffler function 143 



p(t) = 



^ r(i 



(1 + aj) 



(20) 



For a = 1 the Mittag-Leffler function has the standard 
exponential form 



Wt 



exp 



Wt 



(21) 



whereas for < a < 1 it interpolates the initial stretched 
exponential form as 



exp 



Wt" 



r(l + a) 



(22) 



and, however, at the long time, the the initial stretched 
exponential behavior turns over to the power-law behav- 
ior 



Wt° 



1 



r(l + a)M^ 



(23) 



In this study, we have focused the behavior of short time 
limit, and we have written the conditional probability 
(20) for the the simplicity as down 



P (t) = exp 



Wt" 



(24) 



using Eq. (22) form. Actually, this result clearly indi- 
cates that the character of the relaxation has stretched 
exponential i.e. KWW form for the disordered complex 
systems. The jump matrix W in Eq. (24) may be repre- 
sented in terms of collision and unit matrix elements as 



W J-1 



(25) 



where A is the relaxation rate, J is the collision matrix, 
and 1 is a unit matrix. 

For the TV-level the relaxation rate is given by 



A = Nw. 

and the collision matrix J is prescribed by 



(26) 



J = 



/ 1/iV 1/iV 
1/N 1/N 



V 1/A^ l/N 



1/N\ 
1/N 



1/N J 



(27) 



where it; is the jump rate from one value of x to another. 
The jump matrix W is also given by 



W : 



/{l~N)w w 

w {1-N)w 



\ 



w 
w 



\ 



(28) 



(1 - N)w J 



The usefulness of the decomposition in Eq. (25) is 
borne out by the fact that the J matrix has a very simple 
property; it is idempotent, i.e.. 



J2 = J, J3 ^ J, 



(29) 



for any integer fc > 0. This property allows us to im- 
mediately construct the conditional probability P{t) in 
matrix form: 



P{t) = exp A (^J- t 



(30) 



Hence, using direct power series expansion and Eq. (29), 
Eq. (30) can be written follow 



F (t) = exp (-At") 1 - J + Jexp(At") 



(31) 



It is convenient to introduce the more general notation 
for basis vectors in order to make easy of calculations. 
Therefore, let us associate a stochastic state \n) [n = 
1,2,3, ...,7V for multi-level jumping process) instead of 
\x) . Thus, it is possible that Eq. (10) is rewritten down 
form 

{x (0) X [t)) ^ {n\ X \n) {m\ P [t) \n) (m| X |m) 

(32) 



with apriori occupation probability 
1 



Pn{x) 



N' 



n = 1,2,. ..,7V. 



(33) 



Now, it is possible expectation values of the quantities 
which we interested in; 

(n| P {t) H = ^+ (Sr^rn - ^) exp (-A^) (34) 
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where n, m = 1, 2, A''. The matrix of the operator X 
is 

{n\ X \m) — XnSnm, n,m = 1,2, N (35) 

where X„ are the allowed values of the stochastic vari- 
ables. In addition, average of the collision matrix J is 
presented as 



{n\ J \m) 



1 



n,m — 1,2, N 



(36) 



If Eqs. (34) and (35) are inserted in Eq. (10), after a bit 
of algebra, the correlation function can be written down 
form 

{x (0) X (t)) = {xf + - (x)') exp i-Xt") (37) 

where the deterministic quantities are weighted averages 
over the available states of the corresponding variable. 
The average value of {x) and {x^) in the stationary state 
are given by 

N N 



and 



N N 

(x2)=5]p„(n|r2|n) = -^X2 (39) 



respectively. 

The response function (t) in Eq. (4) for this process 
is obtained as 



* (<) = ^ {(^') - {^f [{^') - i^f] exp i-Xn] 

(40) 

using Eq. (37), on the other hand, the limit the behavior 
of '5 (i) at the infinity is presented as 



(t = CX)) 



kT 



(x^) {xY 



(41) 



When Eqs. (40) and (41) are inserted into Eq. (3), the 
relaxation function can be obtained down form 



<^{t)^^{{Ax{tyf)exp[~Xn (42) 



where y(Aa;) y is the mean squared displacement which 
contains solely the molecular contributions 



{Ax{t)f)^{x^)-{xy 



2Kr. 



r(l + a) 



(43) 



where is known the diffusion constants which is a gen- 
eralization of the Einstein- Stokes- Smoluchowski relation 
1^. The Ka is defined as 

Kc, = kBT/mr]^ (44) 
where rja is the friction coefficient which is a measure for 
interaction of the particle with its environment, and ni 
denotes mass of the particle H S IH 113 . 

Finally the relaxation function in Eq. (42) can be sim- 
plified as 



1 



2K^ 



kTr{l + a) 



e exp [-Ar] 



(45) 



The prefactor in Eq. (45) can be taken $o which is a con- 
stant in Eq. (2). As a result, the relaxation function (45) 
is consistent with Eq. (2), which clearly indicates that 
fluctuation quantity in a disordered complex system de- 
cays with time as stretched exponential i.e, KWW form. 

In this study, applying the multi-level jumping formal- 
ism to the fluctuation quantity which make diffusive mo- 
tion stochastically in the disordered complex system, we 
have analytically obtained the relaxation function KWW 
form in terms of correlation function in absence of the ex- 
ternal field. It is concluded that multi-level jumping pro- 
cess formalism is quite powerful technique for the mod- 
elling of the Brownian motion to obtain the relaxation 
function of disordered complex systems. 



[1] J. C. Phillips, Rep. Prog. Phys. 59, 1133 (1996), and 

references there in. 
[2] P. Debye, Polar Molecules (Chemical Catalog, New York, 

1929, reprinted by Dover Publications, New York, 1954). 
[3] R. Kohlrausch, Ann. Phys. (Leipzig) 12, 393 (1847). 
[4] G. Williams and D. C. Watts, Trans. Faraday Soc. 66, 

80 (1970). 

[5] E. W. Montroll, G. H. Weiss, J. Math. Phys. 6 167 (1965). 
[6] H. Scher, M. Lax, Phys. Rev. B 7 4491 (1973). 
[7] E. W. Montroll and H. Sher, J. Stat. Phys. 9 101 (1973). 
[8] G. Weiss and R. Rubin, Adv. Chem. Phys. 52 363 (1983). 
[9] R. Metzler and J. Klafter, Adv. Chem. Phys. 116, 223 
(2001); Phys. Rep. 339, 1 (2000); J. Non-Cryst. Sol. 305 

81 (2002); J. Molec. Liq. 86, 219 (2000). 



[10] J. Klafter, A. Blumen and M. F. Shlesinger, Phys. Rev. 

A 35 3081 (1987). 
[11] J. Klafter and R. Silbey, Phys. Rev. Lett. 44 55 (1980). 
[12] B. Movaghar, W. Schirmacher, J. Phys. C textbfl4 859 

(1981). 

[13] O. Bleibaum, H. Bottger, V.V. Bryksin, Phys. Rev. B 56 
6698 (1997). 

[14] T. Wichmann, K. W. Kehr, J. Phys.: Condens. Matter 

7 717 (1995). 
[15] S. H. Glarum, J. Chem. Phys. 33 639 (1960). 
[16] S. R. Elliott, A. P. Owens, Philos. Mag. B 60 777 (1989). 
[17] J. Klafter, A. Blumen, Chem. Phys. Lett. 119 377 (1985). 
[18] M. F. Shlesinger, E. W. Montroll, Proc. Natl. Acad. Sci. 

USA 81 1280 (1984). 



5 



[19] H. Sher and E. W. MontroU, Phys. Rev. B 12 2455 
(1975). 

[20] P. Grassberger and I. Procaccia, J. Chem. Phys. 77, 6281 
(1982). 

[21] K. Funke, I. Riess, Z. Phys. Chem. Neue Folge 140 217 

(1984) ; K. Funke, Z. Phys. Chem. Neue Folge 154 251 
(1987). 

[22] S. Kirkpatrick, Rev. Mod. Phys. 45 574 (1973). 

[23] P. Soven, Phys. Rev. 156 809 (1967). 

[24] T. Odagaki, M. Lax, Phys. Rev. B 24 5284 (1981). 

[25] I. Webman, Phys. Rev. Lett. 47 1496 (1981). 

[26] R. G. Palmer, D. L. Stein, E. Abrahams and P. W. An- 
derson, Phys. Rev. Lett. 53, 958 (1984). 

[27] D. Kumar and S. R. Shenoy, Phys. Rev. B 34 3547 
(1986). 

[28] A. T. Ogielski and D. L. Stein, Phys. Rev. Lett. 55 1634 

(1985) . 

[29] S. Teitel and E. Domany, Phys. Rev. Lett. 55 2176 
(1985). 

[30] L A. CampbeU, J. Phys (France) Lett. 46, L1159 (1985). 
[31] P. Jund, R. Jullien and LA. CampbeU, Phys. Rev. E 63 

036131 (2001). 
[32] A. Perez-Madrid, Phys. Rev. E 69, 062102 (2004). 
[33] M. H. Vainstein, I. V. L. Costa, R. Morgado and F. A. 

Oliveira, E-print: cond-mat/0501522 
[34] K. L. Ngai, Comments Sohd State Phys. 9 127 (1979); 

Comments Solid State Phys. 141 (1979). 
[35] E. Aydiner, Phys. Rev. E 71, 046103 (2005). 
[36] Y. P. Kalmykov, Phys. Rev. E 70, 051106 (2004). 
[37] W. T. Coffey, Y. P. Kalmykov, and S. V. Titov, Phys. 



Rev. E 67, 061115 (2003). 
[38] Y. P. Kalmykov, W. T. Coffey, D. S. F. Crothers, and S. 

V. Titov, Phys. Rev. E 70, 041103 (2004). 
[39] W. T. Coffey, Y. P. Kalmykov, and S. V. Titov, J. Chem. 

Phys. 116, 6422 (2002). 
[40] P. G. Debenedetti and F. H. Stillinger, Nature (London) 

410, 259 (2001). 
[41] H. Frauenfelder, S. Sligar and P. G. Wolynes, Science 

254, 1598 (1991). 
[42] R. Hilfer, Fractals 3, 211 (1995). 

[43] R. Hilfer and L. Anton, Phys. Rev. E 51, R848 (1995). 
[44] N. G. van Kampen, Stochastic Processes in Physics and 

Chemistry (North-Holland, Amsterdam, 1981). 
[45] B. D. Huges, Random Walks and Random Environments, 

Vol.1 (Clarendon Press, Oxford, 1995). 
[46] S. Dattagupta, Relaxation Phenomena in Condensed 

Matter Physics (Academic Press, Orlando, 1987). 
[47] K. B. Oldham and J. Spanier, The Fractional Calculus 

(Academic Press, New York, 1974). 
[48] R. Metzler, W. G. Glockle and T. F. Nonnenmacher, 

Physica A 211, 13 (1994). 
[49] A. Erdelyi, Tables of Integral Transforms (McGraw-Hill, 

New York, 1954). 
[50] R. Metzler and T. F. Nonnenmacher, Chemical Phys. 

284, 67 (2002). 

[51] R. Metzler, E. Barkai and J. Klafter, Europhys. Lett. 46, 

431 (1999); Phys. Rev. Lett. 82 3563 (1999). 
[52] E. Barkai, Phys. Rev. E 63, 046118 (2001). 



